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Tidal Love numbers of a slowly spinning neutron star 
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By extending our recent framework to describe the tidal deformations of a spinning compact 
object, we compute for the first time the tidal Love numbers of a spinning neutron star to linear 
order in the angular momentum. The spin of the object introduces couplings between electric and 
magnetic distortions and new classes of spin-induced (“rotational”) tidal Love numbers emerge. 

We focus on stationary tidal fields, which induce axisymmetric perturbations. We present the 
perturbation equations for both electric-led and magnetic-led rotational Love numbers for generic 
multipoles and explicitly solve them for various tabulated equations of state and for a tidal field with 
an electric (even parity) and magnetic (odd parity) component with f = 2, 3, 4. For a binary system 
close to the merger, various components of the tidal field become relevant. In this case we find 
that an octupolar magnetic tidal field can significantly modify the mass quadrupole moment of a 
neutron star. Preliminary estimates, assuming a spin parameter x ~ 0.05, show modifications > 10% 
relative to the static case, at an orbital distance of five stellar radii. Furthermore, the rotational 
Love numbers as functions of the moment of inertia are much more seusitive to the equation of 
state than in the static case, where approximate universal relations at the percent level exist. For 
a neutron-star binary approaching the merger, we estimate that the approximate universality of 
the induced mass quadrupole moment deteriorates from 1% in the static case to roughly 6% when 
X ~ 0.05. Our results suggest that spin-tidal couplings can introduce important corrections to the 
gravitational waveforms of spinning neutron-star binaries approaching the merger. 

PACS numbers: 04.20.-q, 04.25.-g, 04.70.Bw, 04.30.-w. 


I. INTRODUCTION 

In a previous paper [T] (henceforth Paper I) we initi¬ 
ated the fascinating study of the tidal deformations of 
a spinning compact object (see also Ref. [2] which ap¬ 
peared at the same time with Paper I). Our framework 
- and that independently developed in Refs. [iia - is 
based on a perturbative expansion in the angular mo¬ 
ment which is valid for slowly spinning objects. When 
applied to the case of spinning black holes, these studies 
have shown that the deformations of the multipole mo¬ 
ments of a Kerr black holes immersed in a tidal field are 
all zercQ In other words, the tidal Love numbers - i.e., 
the deformation of the multipole moments per unit tidal 
field [H [S] - of a spinning black hole are zero, as in the 
static case [MU- 

In this follow-up, we focus on the tidal Love numbers 
of a spinning neutron star (NS). The main motivation 
for such analysis comes from the prospects of measuring 
the tidal Love numbers through gravitational-wave (GW) 
detections of compact binaries ISEQ] (see also [TTH2T] ). 

The tidal Love numbers encode the deformability of 
a self-gravitating object immersed in a tidal environ¬ 
ment and depend sensibly on the object internal struc¬ 
ture H [g. For static NSs, these numbers depend on 
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the NS mass and on the equation of state (EoS) of the 
matter composing the star. Therefore, measuring the NS 
deformability through GW detections would help to con¬ 
strain the behavior of matter at ultranuclear density [su¬ 
ng. 

Previous studies on the tidal deformability of NSs have 
considered only nonspinning objects. The scope of the 
present paper is to quantify the importance of spin cor¬ 
rections to the NS tidal deformability by computing ex¬ 
plicitly a new class of tidal Love numbers that was intro¬ 
duced in Paper I and in Ref. [2]. 


Executive summary 

For the reader’s convenience, we summarize here the 
structure of the paper and our main results. In Sec. [IT] we 
introduce the tidal Love numbers of a rotating body, ex¬ 
tending the formalism developed in Paper 1. Section |ll B 
is devoted to an estimate of spin effects for the tidal Love 
numbers, anticipating the main results obtained in the 
rest of the paper [cf. Fig. [^. Our estimates show that 
- for a NS-NS binary approaching the mereer - the mass 
quadrupole moment of one star with spirjjx ~ 0-05 can 
deviate roughly by 13% relative to the static case [cf. 
Eq. (§]. 


^ Here y = J/Af^ is the dimensionless spin parameter, M and J 
being the mass and angular momentum of the star, respectively. 
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The equations governing the tidal perturbations are 
discussed in Sec. |III| and in the Appendixes and 
In Sec. |III E| and in Appendix we present the full ex¬ 
terior solution of a tidally deformed spinning object to 
first order in the spin for both an electric and a mag¬ 
netic component of the tidal field with £ = 2,3,4. Such 
an exterior solution is specialized to the case of a spin¬ 
ning black hole in Appendix where we also show that 
the corresponding tidal Love numbers are precisely zero 
in this case. This result complements the analyses pre¬ 
sented in Ref. [3] and in Paper I. 

The final expressions for the rotational Love numbers 
of a spinning NS and the matching procedure to explic¬ 
itly compute them are presented in Sec. [IV] Section |V] 
is devoted to a detailed numerical analysis of the rota¬ 
tional Love numbers in various cases and contains the 
main results of this work. In Fig. [^we show some of the 
rotational Love numbers as functions of the compactness 
for various EoSs. In Fig.|^we discuss to which level these 
new Love numbers are independent of the EoS, showing 
that the ^proximate universality [26U28] that exists in 
the static]^ case deteriorates as the spin increases. 

In the main text we focus on those rotational Love 
numbers which are associated with perturbations that do 
not break the reflection symmetry on the equatorial plane 
(and are therefore more relevant phenomenologically). 
The case of equatorial-symmetry breaking is briefly dis¬ 
cussed in Appendix Along the way we clarify some 
questions related to the tidal Love numbers, for exam¬ 
ple in Appendix we reveal some inconsistency in the 
computation of the magnetic tidal Love numbers that 
appeared in the literature. Through this work, we use 
geometrized G = c = 1 units. Spacetime indices are de¬ 
noted by Greek letters, while space indices are denoted 
by Latin letters, i.e., is the time coor¬ 

dinate associated with the unperturbed compact object. 


II. TIDAL LOVE NUMBERS OF A SPINNING 
OBJECT 


A. Definition 

Let us consider a (generically spinning) compact ob¬ 
ject immersed in a generic tidal environment. We adopt 
the same decomposition of the tidal field as in Ref. [B]. 
Namely, we define the symmetric and trace-free elec¬ 
tric and magnetic tidal multipole moments as Ea^ ai_ = 
[(£- 2)!]-l(GoaiOa,;a3...a,) and6„,...a, = [§(£+!)(£- 
2)!]“^(eaihcG^^0;a3...aj): where Cabcd is the Weyl tensor. 


a semicolon denotes a covariant derivative, Cabc is the 
permutation symbol, the angular brackets denote sym- 
metrization of the indices and all traces are removed. 
It can be shown that the electric (respectively, magnetic) 
moments Eai...ai (respectively, Bai...ae) can be decom¬ 
posed in a basis of even (respectively, odd) parity spher¬ 
ical harmonics with the symmetry axis aligned with the 
spin vector |31 [B] . We denote by Em and Bm the am¬ 
plitude of the electric and magnetic components of the 
external tidal field with harmonic indices (£, m), where m 
is the azimuthal number. The structure of the external 
tidal field is entirely encoded in the coefficients Em and 
Bm (cf. Ref. [6] for details). 

As a result of the external perturbation, the mass and 
current multipole momenttj^(M£ and Sg, respectively) of 
the compact object will be deformed. In linear perturba¬ 
tion theory, these deformations are proportional to the 
applied tidal field. In the nonrotating case mass (cur¬ 
rent) multipoles have even (odd) parity, and therefore 
they only depend on electric (magnetic) components of 
the tidal field; we can define the standard tidal Love num¬ 
bers as 
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and these quantities are independent of the tidal field. 
The definitions above agree, modulo some normalization 
factor, with the standard ones niziiin], usually denoted 
as electric and magnetic tidal Love numbers, respec¬ 
tively. Because the object is spherically symmetric, the 
azimuthal number m is degenerate. Furthermore, parity 
and the angular momentum number i are conserved: an 
electric (i.e. even parity) tidal field with harmonic index 
£ can only deform the mass multipole moment of order 
£, whereas a magnetic (i.e. odd parity) tidal field with 
harmonic index £ can only deform the current multipole 
moment of order £. 

Such degeneracy is broken when the central object is 
spinning. In such case there exist some selection rules 
which are discussed in detail in Paper I and are schemat¬ 
ically depicted in Fig. Each arrow in this figure de¬ 
notes a coupling of linear order in the spin. Therefore, 
to first order in the spin, electric (respectively, magnetic) 
tidal perturbations with multipolar index £ source mag¬ 
netic (respectively, electric) deformations with multipo¬ 
lar index £ ± 1 and vice versa. In particular, the mass 
quadrupole (£ = 2, electric) moment of the star acquires 
0(x) corrections proportional to the magnetic octupolar 
(£ = 3) tidal field. 

Furthermore, a nonaxisymmetric tidal field with index 
£ would deform the multipole moment of the same parity 


® It is important to stress that the tidal Love numbers entering 
the nearly universal relations found in Refs. 12644281 are those of 
nonspinning NSs, whereas the moment of inertia and the mass 
quadrupole moment refer to rotating stellar configurations. The 
spin corrections to the tidal Love numbers of a NS are computed 
in this work for the first time. 


^ We consider the Geroch-Hansen multipole moments 1291 1,301 . 
which are equivalent m to the multipole moments defined by 
Thorne m using asymptotically mass-centered Cartesian coor¬ 
dinates. 
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the one given in Paper 1. Such generalization is necessary 
to accommodate the presence of a magnetic tidal field or 
of a tidal component with i > 2; both possibilities were 
not considered in Paper 1. 

Because of the spin selection rules mentioned above 
(see also Paper I and Ref. [33]), the Love numbers in 
Eqs. and ^ enjoy various interesting properties (cf. 
Fig.[^ 

1. At zeroth order in the spin, the Love numbers in 
Eq. ([^ vanish, whereas the Love numbers in Eq. ([^ 
are different from zero only when (.' = £ and reduce 
to those defined in Eq. Q. As mentioned above, 
the azimuthal number m is degenerate in this case. 


FIG. 1. (color online). Scheme of the spin-tidal coupling in 
the slow-rotation approximation to second order in the spin 
for the electric-led system (top diagram) and for the magnetic- 
led system (bottom diagram). Vt and At generically denote 
polar (i.e., even parity, or electric) and axial (i.e., odd parity, 
or magnetic) perturbations with harmonic index i. The quan¬ 
tity enclosed in a box denotes the component of the external 
tidal field. Perturbations with different parity and harmonic 
index are coupled to the external perturbations by arrows, 
each arrow (taken in either direction) denoting a coupling of 
linear order in the spin. For example, in the top diagram Vi 
sources deformations to second order in the spin, since 
two arrows are needed to go from Vt to 'Pi+ 2 . Purple hori¬ 
zontal arrows denote Zeeman-like couplings which are nonzero 
only in the nonaxisymmetric case, m 7 ^ 0 (cf. Paper I and 
the review [33] for details.) 


and with the same harmonic index £ to first order in the 
spin through a Zeeman-like splitting term proportional to 
the azimuthal number m and denoted with a horizontal 
purple line in Fig. [^ 

Finally, to second order in the spin, an electric (re¬ 
spectively, magnetic) tidal field with multipolar index 
£ sources electric (respectively, magnetic) deformations 
with multipolar indices £ and £ ± 2. Thus, for instance, 
the mass quadrupole moment acquires a O(x^) correction 
proportional to the electric quadrupolar tidal field. 

On the light of such selection rules, we can define a 
more generic set of tidal Love numbers as follows: 
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The Love numbers with the “-I-” subscript measure how 
a mass (respectively, current) multipole moment is de¬ 
formed by the presence of an electric (respectively, mag¬ 
netic) tidal field, i.e. by a field of the same parity. On 
the other hand, the Love numbers with the ” sub¬ 
script measure how a mass (respectively, current) multi¬ 
pole moment is deformed by the presence of a magnetic 
(respectively, electric) tidal field, i.e. by a field of the op¬ 
posite parity. The definition just introduced generalizes 


2. At first order in the spin, the Love numbers in 
Eq. © are different from zero only when £' = £ ± 1, 
whereas the Love numbers in Eq. ([^ are different 
from zero only when £' = £ and in the nonaxisym¬ 
metric case, m 7^ 0. These latter corrections corre¬ 
spond to the Zeeman splitting mentioned above. 

3. At second order in the spin, the Love numbers in 
Eq. (§ acquire new corrections only in the non¬ 
axisymmetric case, whereas the Love numbers in 
Eq. (§ are different from zero both when £' = £ 
and when £' = £ ± 2 for any value of m. The terms 

and include O(x^) corrections to the 

static Love numbers in Eq. Q. 

In this work we focus on axisymmetric tidal perturba¬ 
tions to first order in the spin. In this case, the Love 
numbers in Eq. (§ reduce to the static case, Eq. 0’ 
and one is left with the spin-induced (“rotational”) Love 
numbers in Eq. ([^ with to = 0. An explicit computation 
of the second-order and nonaxisymmetric corrections is 
left for future work. To simplify the notation, in the 
following we will introduce the shorthand notation 

^ xf ';^, 5X^1!'^ ^ x%l ^. ( 4 ) 

For example, denotes the correction to A^^ arising 

from the coupling to the axisymmetric magnetic octupo- 
lar (£ = 3) component of the tidal field to first order in 
the spin. Among the various rotational Love numbers 

( 23 ) 

that we will compute, i5A^ ' is the most relevant one 
from a phenomenological point of view because it is as¬ 
sociated with a deformation of the quadrupole moment 
of the star. 


B. Estimating spin corrections to the tidal Love 
numbers 

Although the distribution of NS spin is uncertain, old 
NSs in the late stages of a binary inspiral are expected to 
rotate rather slowly. The fastest spinning NS observed 
so far is the most massive companion of the double pul¬ 
sar system PSR J0737-3039A [31], with a spin period 
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FIG. 2. (color online). Ratio between the (dimensionless) 
rotational tidal Love number and the standard electric 

quadrupolar tidal Love number for the various EoSs con¬ 
sidered in this work. This ratio is crucial for the estimate ([^. 


of « 23 ms, which corresponds to x 0.02 — 0.05, de¬ 
pending on the EoS [35113S] ■ Such rotation rate is not 
expected to decrease substantially as this system ap¬ 
proaches the merger (see Ref. m for a recent discussion). 
Furthermore, the observation of numerous isolated mil¬ 
lisecond pulsars suggests that higher spin rates (as high 
asx ~ 0.1 [37]) might alos be found in NS-NS binary sys¬ 
tems. In light of this uncertainty, our estimates in this 
paper are based on the assumption x ~ 0.0^ 

It is therefore relevant to discuss in which regime the 
0{x) spin corrections computed in this paper are phe¬ 
nomenologically relevant. To answer this question, let us 
consider - as a representative example - an equal-mass 
NS-NS binary with total mass ^ 2M at large orbital 
separation tq ^ M. To the leading order, the £ = 2 
components of the tidal field scale as (cf. e.g. Ref. [38] ) 


£( 2 ) _ — 
^rn 3 

'o 



(5) 


(modulo a prefactor that depends on m) for the electric 
component and the magnetic component, respectively. 
Note that the la tter component is suppressed by a fac¬ 
tor V = \/MItq relative to the electric component. On 
the other hand, the I = "£> components involve a covari¬ 
ant derivative of the tidal tensor [S] and are therefore 
suppressed by a factor l/rp relative to their t = 2 coun¬ 
terparts, i.e. 


£(3) ^ — 

4 ) 

'o 


H V Tq 


( 6 ) 


As discussed above, spin corrections introduce new 
couplings between the multipole moments of a spinning 


NS and various components of the external tidal field. 
In Sec. IV A , we will show that the deformed quadrupole 
moment of a spinning NS in an axisymmetric tidal envi¬ 
ronment to linear order in the spin reads 



where tilded quantities are made dimensionless and 
spin independent by dividing the corresponding quan¬ 
tity by suitable powers of M and x [as explained below 
Eqs. ([3^ -(38)]. As we will show, relations similar to 
Eq. Q can be obtained also for other (mass and current) 
mnltipole moments and involve the novel family of tidal 
Love nnmbers defined in Eqs. (§-(§. The main pur¬ 
pose of the present work is to compute such spin-induced 
corrections for a NS described by a realistic EoS to first 
order in the spin. 

To quantify the effect of such spin-induced corrections, 
we can estimate the 0(x) axisymmetric tidal correction 
to M 2 relative to the static case. Namely, from Eq. 0 
we obtain 


5M2 


13% 


5Xf^\ 



where is the induced quadrupole moment Q for 

X = 0, 5M 2 is its 0{x) correction, C = M/R is the stellar 
compactness, and we have normalized the parameters to 
X ~ 0.05, C ~ 0.2 and rg ~ 5R. The latter represents an 
extrapolation because the expression above is only valid 
when the source of the tidal field is at large distance from 
the central object, rg » R, but it should nevertheless 
capture the correct order of magnitude. 

The estimate Q is one of the main results of this 
paper. It is clear that the relative amplitude depends 
on the ratio which we show in Fig. ^ as a 

function of the NS compactness, anticipating the main 
computation presented in the rest of this work. Inter¬ 
estingly, for a typical NS compactness, M/R ks 0.2, 
JA^^^ « 30A^^ and, from Eq. ([^, we find that, for 
X « 0.05, the 0{x) correction to M 2 is about 13% of 
its static counterparl]^ Because the amplitude of higher- 
order tidal multipoles increases near the merger, we ex¬ 
pect such correction to become even larger. This simple 
estimate shows that spin corrections to the tidal Love 
numbers are non-negligible and they might considerably 
affect the gravitational waveforms near the merger. 

Finally, another relevant question concerns second- 
order spin corrections. The relative amplitude of first- 


® Obviously, because we consider linear perturbations in the spin, 
the magnitude of the corrections discussed in the rest of the paper 
would be double in the unlikely case of a spinning NS in a binary 
system with x ~ OT near the merger. 


Our estimates take into account only the axisymmetric compo¬ 
nents of the tidal field. Since the background is axisymmetric, 
perturbations with different azimuthal number m are decoupled 
and nonaxisymmetric perturbations would simply introduce ex¬ 
tra terms in Eq. 0 which are proportional to new Love numbers. 
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and second-order spin corrections to the tidal Love num¬ 
bers depends on two competitive effects: on the one hand 
the corrections linear in the spin are larger than the 
quadratic corrections when x ^ 1, as expected for old 
NSs in the late stages of a binary inspiral; on the other 
hand, the magnetic octupolar component of the tidal field 
is suppressed relative to the quadrupolar electric compo¬ 
nent, at least when the orbital velocity is small, n <C 1. It 
is therefore possible that 0{x^) spin corrections can also 
be relevant even in the x <C 1 regime, at least for inter¬ 
mediate orbital distances. Clearly, in order to quantify 
the relative magnitude of second-order spin corrections 
(as well as of nonaxisymmetric 0(x) corrections), it is 
necessary to compute other rotational Love numbers (for 
instance ). We leave this important task for future 

work, focusing here on first-order, axisymmetric correc¬ 
tions. 


III. PERTURBATION EQUATIONS 


where P = P{r) and p = p{r) are the background pres¬ 
sure and energy density, respectively. To hrst order in the 
spin, = e“'^/^(l, 0, 0, Q) where fl is the fluid angular 
velocity and the normalization constant ensures = — 1. 
Einstein’s equation^ = 2kT^„, reduce to the classi¬ 
cal Tolman-Oppenheimer-Volkoff equations for the back¬ 
ground variables 


, 2(M + Kr^P) 


P' = - 


r{r — 2Ai) ’ 

{P + p) (At -I- Kr^P) 
r(r — 2A4) 


M! = nr'^p , 


Kr{P + p) {ruj'+ Aoj) 4_, 


r-2M 


- OJ 

r 


( 11 ) 

( 12 ) 

(13) 


where we dehned w = Q — w. Assuming a barotropic 
equation of state in the form P = P{p), these equations 
can be integrated numerically with standard methods. 


The study of linear perturbations of a slowly rotating 
star has been initiated in Refs. and was recently 

extended in the context of black-hole perturbation the¬ 
ory (see [35] for a review). Here we follow exten¬ 

sively the formalism developed in Paper I, to which we 
refer for further details. In particular, we make the same 
working assumptions discussed in Paper I; namely, we 
consider a slowly rotating object, whose spin axis coin¬ 
cides with the axis of symmetry of the tidal field, which is 
assumed to be axisymmetric. The tidal field is weak and 
varies slowly in time, and in practice we consider only 
stationary tidal perturbations. As we will show, axisym- 
metry is in fact required for consistency of the held equa¬ 
tions in the stationary case. As explained in Paper I, our 
approximations are valid for a compact binary at a large 
orbital distance. Furthermore, NSs in a binary system 
are expected to rotate slowly so that the slow-rotation 
approximation is well justihed (see, e.g., [361137]). At 
variance with Paper I, we limit our analysis to hrst order 
in the spin, but consider both an electric and a magnetic 
tidal component with generic multipole 


A. Spinning background 


We consider a slowly rotating object whose line ele¬ 
ment, to hrst order in the spin, can be written as mi 03 


ds^= —e'^dt^ + 


dr'^ 


l-2Mlr 


— 2 sin^ 'duir'^dtdip + r^dVL^ 


(9) 


where = d{P + sin^ ddip^, and the radial functions v 
and M are of zeroth order in rotation, whereas ui is of 
hrst order. The stress-energy tensor of the perfect huid 
is 


( 10 ) 


B. Tidal perturbations 


Metric and huid perturbations are decomposed in 
spherical harmonics according to their parity as dis¬ 
cussed in Appendix [^ Since we are considering sta¬ 
tionary tidal helds, we assume no explicit time depen¬ 
dence in the perturbations and no huid motiorj^ setting 
= 0 in the huid perturbations 6u^, 

cf. Eqs. (|^-(|A^. 


By using this decomposition, solving the gravitational 
equations perturbatively, and making use of the orthog¬ 
onality properties of the spherical harmonics [31 iO], 
we obtain a system of ordinary differential equations for 
the radial variables where various spin couplings between 
diherent multipolar indices £ and between perturbations 
with opposite parity are present. 

The equations separate into two groups according to 
their parity. In Paper I, we dehned the so-called polar-led 
and axial-led systems, which describe the metric defor¬ 
mations induced by a purely electric and a purely mag¬ 
netic tidal held, respectively. It can be shown that con¬ 
sistency of the held equations requires stationary pertur¬ 
bations to also be axisymmetric. In other words, the 
absence of explicit time dependence in the perturbations 
given in Appendix [^ also imposes m = 0, where m is 


^ In our geometrized units k = 47r but we will leave it unspecified 
to keep track of the coupling to matter in the equations. 

® Recently, it has been noted that an irrotational (rather than a 
static) fluid provides a more realistic configuration for a station¬ 
ary tidally distorted object in a binary system SSj. Such condi¬ 
tion only affects the odd-parity sector of the perturbations, but 
has dramatic consequences for the magnetic tidal Love numbers. 
Here we consider only static fluids, this assumption is consistent 
with the field equations to first order in the spin. An extension 
of the work of Ref. m to account for spin corrections is left for 
the future. 


= (P + p)uf^u'^ + , 
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the azimuthal number of the perturbations, which is con¬ 
served by virtue of the axisymmetry of the background. 
On the contrary, tidal perturbations with m 7 ^ 0 must 
have a weak time dependence (cf. e.g. Ref. [3] and Pa¬ 
per I). Therefore, to first order in the spin and in the 
stationary case, the polar-led and axial-led systems given 
in Paper I schematically reduce to 


Pl = 0 


Al+ 1 + iaQL+lPL = 0 , 

(14) 

Al- 1 + ^aQhVL = 0 


Al = 0 

'Pl+I + ^aQL+lAL = 0 , 
Vl-I + CaSidi = 0 

(15) 


where Cq is a bookkeeping parameter for the order of 
the spin perturbations, Ai^ Ai (respectively, Vt, Pi) are 
linear combinations of the axial (respectively, polar) per¬ 
turbations with multipolar index ^ and Qi = £/\/AP — 1 
are related to the Clebsch-Gordan coefficients in quan¬ 
tum mechanic^ The first set of equations in Eqs. (14) 
and © is of zeroth order in the spin and corresponds 
to the well-studied static case whereas the 

remaining two sets of equations in Eqs. (14) and (151 are 
spin-induced corrections. As clear from the structure of 
the above equations, polar perturbations with multipolar 
index t source axial perturbations with multipolar index 
£ -|- 1 and £ — 1 to first order in the spin and vice versa, 
cf. Paper I for details. 

Remarkably, it is possible to reduce the systems above 
to a simple set of differential equations for a generic mul¬ 
tipole £. To simplify the notation, we will expand any 
variable to first order in the spin, e.g. and 

will denote a function of zeroth order in the spin and its 
first-order correction, respectively. 


C. Electric-led tidal perturbations 


For any £ > 2, the polar-led system (141 reduces to the 


following coupled ordinary differential equations (ODEs) 


V, 

V' 

V 




H, 


W 


(t-ti) 

* 

(t-i) 


= 0 

G+1)' 

0 

{£-!)' 


_ g*. (£) 


(16) 


where we defined the differential operatorj^ 


V 

V 




(^) 


I 


(«) 




(17) 

(18) 


The coefficients and {i = 1,2) and the source 

terms are given in Appendix B As expected, the 

(t) . ' - ' 

source terms depend on Hq , i.e. polar perturbations 
with multipolar index £ source axial perturbations with 
£ ± 1. The remaining variables are algebraically related 
to and by relations given in Appendix [b| 

In particular, = 0 = so that the perturbed 

metric remains invariant under t ^ —t and if ^ —tp. 
The perturbatio ns also depend on the speed of sound in 
the fluid, Cs = y/dPjdp. 

Note that, in principle, the equations governing axial 
perturbations with £ = 1 are different from those with 
£ > 2. When £ = 1, one can use a residual gauge freedom 
to set = 0 = Sh^'^ [^[^. In the stationary case, it 
turns out that = 0 for any £ so that the equations 

above also include the special case £ = 1. As we discuss 
below, this property does not hold in the magnetic-led 
case. 


D. Magnetic-led tidal perturbations 


For any £ > 2, the axial-led system (15) reduces to the 
following coupled ODEs 


( 

So 

= 0 

1 


[ %-i) 

1 1 

1 




= 5'. 


(A 


(19) 


where again the source terms are given in Ap¬ 

pendix]^ and depend on and its derivatives. In this 
case, axial perturbations with multipolar index £ source 
polar perturbations with £ ± 1 . 

Polar equations with £ = 1 are governed by a different 
set of equations and must therefore be discussed sepa¬ 
rately. When £ = 1, one can set = 0 = Sh[P in 
the metric ansatz without loss of generality. For this rea¬ 
son, the differential operator and the source term 
5 ( 1 ) 

are different from those presented above and their 
explicit form is given in Appendix [B| 


^ In the context of quantum mechanics, the spin selection rules are 
enforced by the addition rules for angular momenta m- 


Henceforth, we will use an asterisk to denote quantities related 
to magnetic perturbations. 
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E. Exterior solution 


Although the set of equations just presented is valid for 
any £ > 2 , finding the vacuum solution without specifying 
the value of £ is rather challenging. Indeed, to zeroth 
order in the spin the polar and axial solutions of the first 
equation in (16) and (19), respectively, can be written in 
terms of special functions [siiTiiini 


= A P!{r/M -1) + B Q\(r/M - 1), 
2 


ABG 


2,0 

2,2 



1 -f,^ + 2 \ 
- 1,2 ) ■ 


( 20 ) 



where Pj and are associated Legendre functions, iF\ 
is the hypergeometric function, 6 * 2,2 is the Meijer func¬ 
tion, and A, B, C, D are integration constants. When 
plugged into the source terms, the corresponding first- 
order equations do not appear to have an analytical so¬ 
lution for generic £. However, the special functions above 
simplify enormously when £ is an integer (i.e. in the case 
of interest) and therefore the exterior perturbation equa¬ 
tions can be solved analytically for any specific integer 
value, ^ = 2,3,4,... 


The explicit vacuum solution to the £ = 2 electric-led 
sector reads (cf. Paper I) 


H, 


= a2y{y-2) - 72 -3-f ^ -f 3j/-f ^{y-2)y\og [1 - 2/y\ 


^fe( 3 )^_ Mx 
° 6720y2 


-128^352/(5?/ - 4)02 + 3 (-16 - 44?/ -h 90?/2 -f 21Qy'^ - 945?/^ -f 405?/'^ 

+y (64 - 80?/ -f 1080?/3 - 1350?/^ -f 405?/®) tanh’ 


1 -y 


72 


+ 35 ?/ (8 + 20?/ + 60?/2 - 210?/® + 90?/^ + 15?/® (8 - 10?/ + 3?/2) log [1 - 2/y\) 732 )] 


= - 


Mx 

30?/2 


- 20 yi 5 ?/®a 2 + yi5 (4 - 30?/" + 30?/® + 3?/ (4 - 10?/" + 5?/®) log [1 - 2/y]) 72 - 30?/7i2 


( 22 ) 


(23) 

(24) 


where y = r/M and 02 , 72 , 723 , 721 are integration con¬ 
stants^ The remaining nonvanishing perturbations can 
be obtained from the solution above using the algebraic 
relations presented in Appendix [B| 

Note that, in the solution above, we have already fixed 
the constants 023 and 021 to avoid spurious components 


of the tidal field which arise from the homogeneous sys¬ 
tem associated with Eq. ( [I^ (cf. Paper I for details). 
With this choice, the external tidal field has a purely 
electric, quadrupolar component. 

Likewise, the vacuum solution to the £ = 2 magnetic- 
led sector reads 


M 


^0 = -^ [-24(y - 2)?/®a; -k (-4 - 4?/ - 6 ?/" -f 6 ?/® -f 3(?/ - 2)?/® log [1 - 2/y]) 7 J] , 


(25) 


= 


__ 

0 “ 120(2/ - 2)2y3 L 


X 


-72v^(2/ - 2)y^a*2 + (2 {A + 2y- 6 ?/" - 23//® - 6ly‘^ + 455//® - 420//® -f I05y'^) 


140(2/ - 2)?/4 L 

+3y* (-6 - 1372 / + 280//" - 175?/® -f 35//^) log [1 - 2/y\) 7 ^ -k 70//® (2 (-2 - IQy + 65//" - 60//® -f l£>y'^) 
+15(2/ - 2 )^( 2 / - 1 ) 2 /' log [1 - 2/y]) 7 * 2 ] , (26) 


-24\/T5(y - 2)"y®(102/ - 9)03 + Vl5 (I 6 - 42/ - 32//" - 2My'^ + 60//® 


+3(2/ - 2 )" 2 /®( 102 / - 9) log [1 - 2/2/]) 72 - 1202 /® 7 i* 2 ] , 


(27) 


where 
As m 




in a%y 7^,1 75 q, 7 q-i, are integration constants, 
^eciric-red dase, Vae other constants aJs and 


a 2 i have been fixed by requiring that the tidal field has 
a purely magnetic, quadrupolar component. 


Our notation is the following. The constants ai and 7 ^ are 
proportional to O(x^) perturbations, whereas are associated 
with 0{x) perturbations with i = i which are sourced hy i = j 


tidal perturbations with opposite parity. As before, we denote 
quantities related to the magnetic sector with an asterisk. 
























The equations above represent the generic exterior so¬ 
lution of a slowly spinning object immersed in an ax- 
isymmetric, quadrupolar tidal field, whose electric and 
magnetic components are proportional to 02 and a^, re¬ 
spectively. This solution was computed in Ref. [2] using 
a different notation and in the generic, nonaxisymmetric 
case. The explicit vacuum solutions for the cases with 
£ = 3 and £ = 4 are instead new and given in Appendix [C| 
To identify the electric and magnetic components of 
the tidal held we can expand the exterior solution at large 
distances. In the axisymmetric case, we can extract £q ' 
and from the asymptotic behavior of the gtt and gt^ 
components, namely [ 6 ] 


gauge-invariant way by evaluating the binding energy of 
a point particle in a circular geodesic of an axisymmet¬ 
ric, asymptotically-flat spacetime [55H55] . The multipole 
moments can then be identihed by comparing the binding 
energy order by order in a (gauge-invariant) low-velocity 
expansion. 


A. Multipole moments of a tidally distorted NS 

Further details on this procedure are given in Paper I, 
here we simply present the hnal result for the hrst mul¬ 
tipole momenta^ 



(28) 


(29) 


where Sf„ is one of the vector spherical harmonics dehned 
in Eq. (A4). The normalization of the dimensionless con¬ 
stants ai and aj in the analytical solution just presented 
has been chosen such that 


M 2 

272 


[ 7573* 

272*3 \ 

\ 

M3 


X 



[ 61 ) 

^3 

373* 





M4 

4^77? 





6x 

' V 5 

-^72 + 74 + 

^ (732 + 734) 

, (32) 

II 

474 

7 ^ 


(177773* - 67274*3) 
11760? 

(33) 


c(£) _ f-{f- - 1 ) Oit 

— 




- 2 M^- 


IV. ROTATIONAL LOVE NUMBERS OE A 
SLOWLY ROTATING NS 


With the explicit exterior solution at hand, we can now 
compute the spin-induced (rotational) Love numbers. As 
explained in detail in Paper I, the constants and a* 
are associated with the i = i electric and magnetic com¬ 
ponents of the external tidal held [cf. Eq. (30)], whereas 
the constants 7 ^ and jij (respectively, 7 * and 7 ,*-) are as¬ 
sociated with the deformation of the mass (respectively, 
current) £ = i multipole momentj^ 

To show this fact explicitly, we adopt the procedure 
explained in Paper I to compute the Geroch-Hansen mul¬ 
tipole moments [liEO] of the tidally distorted object. In 
brief, the procedure is based on considering the asymp¬ 
totically fiat solution of the metric perturbations, set¬ 
ting = a* = 0 in the equations above. Using this 
solution, the multipole moments can be computed in a 


The separability between external tidal field and object’s re¬ 
sponse suffers from some subtle ambiguity that we discussed in 
detail in Paper 1. Here we will follow the prescription of that pa¬ 
per, which seems robust enough at least to first order in the spin. 
In practice, because the Love numbers scale with high powers of 
the inverse compactness, ~ (R/M) « 001 IO], where the power q 
depends on i and on the parity (for example ~ (R/M)^, see 
Fig. El below) such ambiguity can at most affect the Love num¬ 
bers by a small offset, ~ {Rbu/^Bh)^ "C (^ns/^Ns)^? with 
little observational consequences. 


With the multipole moments at hand, the Love num¬ 
bers in Eqs. Q and © can easily be computed. By using 
Eq. (30), to zeroth order in the spin we obtain 


(2) _ 272 (3) _ 73 ^-^’’ ,.( 4 ) _ 

^ V^a2’ " ^ 

(34) 


whereas the new 0{x) rotational Love numbers read 


6 (7573573*-22472*3) 

(35) 

^ 28777ra2 

(36) 

5^(34) .^83(4474 + 77734) 

(37) 

eA43) . 3(177773* -67274*3) 

' 105847+a* ■ 

(38) 


The dimensionless quantities previously introduced are 
defined by dividing the expressions above by suitable 
powers of x and M" (where, e.g., n = 6 for ' and 
SX^^\ whereas n = 8 for and ^A^^^) in order 

to make them pure numbers. For example, (5A^ = 


The tidal perturbations also source a deformation of Si, which 
can be reabsorbed in the definition of the NS angular momentum. 
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In terms of these dimensionless quantities, the axisym- 
metric deformation of the mass quadrupole moment to 
linear order in the spin can be expressed as 


Ml - x(2)c(2) 
M3 " 




(39) 


We stress that the above expression contains only the 
axisymmetric components of the tidal perturbations; the 
full expression to first order in the spin would also con¬ 
tain extra terms associated with m ^ 0 perturbations. 
Clearly, to estimate Eq. (39), we need to evaluate the ax¬ 
isymmetric components of the quadrupolar electric tidal 
field and of the octupolar magnetic tidal held. For con¬ 
creteness, let us consider an equal-mass NS-NS binary 
with total mass ~ 2M in a quasicircular motion at or¬ 
bital distance tq. The electric tidal tensor Sab and the 
magnetic tidal tensor Babe in such conhguration were 
computed in Ref. [56]. By comparing the large-distance 
behavior of the perturbed black-hole metric derived in 
Ref. (cf. Eqs. (3.2)-(3.3) in that paper) to our 
leading-order expressions (28)-(|^, we obtain, to New¬ 
tonian ordeJ^'^ 


( 2 ) _ 

^0 - „3 

'o 


( 3 ) _ 18M 

"o - ^4 


M\ 

roj 


1/2 


(40) 

which complement Eqs. © and by including the cor¬ 
rect prefactors for a circular equal-mass binar}]^ Fi¬ 
nally, by inserting Eq. ( |40| into Eq. (39), we obtain the 
expression anticipated in Eq. 0 - Xs previously dis¬ 
cussed, by evaluating Eq. ([^ and using the numerical 
results presented in Sec. [V] below, we estimate that the 
axisymmetric 0{x) spin correction to the quadrupole mo¬ 
ment M 2 for X « 0.05 is about 13% when tq ~ 5R, cf. 
Eq. (|^. 


So far we have focused on perturbations which do not 
break the rehection symmetry with respect to the equato¬ 
rial plane at i? = tt/2. However, an axisymmetric electric 
(respectively, magnetic) tidal held with odd (respectively, 
even) values of i generates perturbations which break the 
equatorial symmetry of the background Such compo¬ 
nents of the tidal held would induce multipole moments 
such as M 3 , S 2 , etc... Our formalism can directly ac¬ 
commodate these new Love numbers, although they are 
only relevant for tidal sources that break the equatorial 
symmetry. 

A binary system of two nonspinning compact objects 
or of two spinning compact objects with (anti)aligned 


Note that, to leading Newtonian order, the perturbed black-hole 
metric computed in Ref. IMI coincides with that of a perturbed 
NS at large distance. For this reason the results of Ref. m are 
directly applicable to our case. 

We note that our expression for coincides with that derived, 
e.g., in Ref. after taking into account their different normal¬ 
ization of the spherical harmonics relative to ours (we use the 
standard definition of [ 57 ]). 


spins enjoys an equatorial symmetry relative to the or¬ 
bital plane. However, in more general conhgurations 
- e.g. if the angular momenta of the two objects are 
not (anti) aligned - the equatorial symmetry can be bro¬ 
ken. This misalignment might happen after the NS 
formation due to supernova kicks, although subsequent 
tidal dissipation can contribute to realign the spins (cf. 
e.g. Ref. [SS] and references therein). For completeness, 
we discuss this problem in Appendix [^ where we also 
present some results for the rotational Love numbers with 
equatorial-symmetry breaking. 


B. Matching and nnmerical procednre 


The above Love numbers can be directly computed 
once the constants 7 ’s are known. The values of these 
constants are fixed by a continuous matching of the nu¬ 
merical interior solution describing the deformed com¬ 
pact star with the analytical exterior solution. Since we 
include only linear corrections in the tidal field, the ra¬ 
tios of the 7 ’s to the corresponding a’s are pure numbers 
independent of the tidal held itself. 

For example, to zeroth order in spin, for the correction 
to the quadrupole mass moment one obtains m 


HZ) _ 
Ae — 


X 


^(1 - 2Cf [{2C - l)p2P -20 + 2] 

V Stt 

[3(1 - 2Cf ((20 - 1 )/ 32 p -20 + 2) log(l - 20) 
+20 ((40-^ + 60^ - 220^ + 150 - 3) /32P 
+40^ - 40^ + 260^ -240 + 6 )] , (41) 


where O = M/R, /3p^ = RH^ evaluated at the 

radius and we recall that tilded quantities are made di¬ 
mensionless by dividing the corresponding quantity by 
suitable powers of M, e.g. = X^e'^/M^. In the small 

compactness limit, one recovers the Newtonian result, 
X^E ~ [ISI- To compute the Love numbers, it is 
therefore sufficient to evaluate the metric perturbations 
at the radius of the object after a numerical integration 
from the interior. The same procedure can be straightfor¬ 
wardly applied to any of the tidal Love numbers dehned 
in Eqs. ([^. 

The case of the rotational Love numbers dehned in 
Eq. ([^ is similar in spirit but slightly more involved. In 
this case the corresponding perturbation equations are 
inhomogeneous ODEs. To match the interior solution 
with the one in the exterior, we adopt the following pro¬ 
cedure. First, we construct a particular solution of the 
inhomogeneous equations which is regular at the origin. 
Then, we also construct the solution to the associate ho¬ 
mogeneous system which is regular at the origin. The 
amplitude of this solution is a free parameter that we 
adjust to match the full interior solution (i.e. a linear 
combination of the inhomogeneous and of the homoge¬ 
neous solutions) to the exterior. It is easy to show that 
imposing continuity at r = i? fixes uniquely the constants 
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7 ’s and also the amplitude of the homogeneous solutions 
for different values of i. In order to increase accuracy, we 
compute higher-order series expansions near the center 
of the star and use such analytical solutions to start the 
numerical integration outwards. 

As a representative example, we present the schematic 
form of the 0 {x) correction to the mass quadrupole mo¬ 
ment: 




Ai + A 2 log(l - 2 c) 

A 3 + Ai log(l - 2c) -I- As log(l - 2c)2 ’ 


(42) 


where Ai are cumbersome polynomial terms that de¬ 


pend on C, //iq 


(3) q{‘^) _ /A 

’ Pph — 


Ph — J 

(' 2 ) 

with being the homogeneous solution associated 

with the third equation of the system (191, and on 

s = RSHq^^ This last term represents the rel- 

f2') 

ative amplitude of the polar perturbation SHq relative 
to its axial source to first order in the spin. All 
these perturbations are evaluated at the radius of the 
star. Note that the value of s might also depend on the 
compactness so that, to explore the Newtonian regime, 
one cannot just take the C —)■ 0 limit of Eq. (|42|. In 


Sec.0 


we will give numerical evidence that ~ l/C® 

in this limit, so that the ratio A^ finite as C —)■ 0 . 

The expressions for the other rotational Love numbers 
defined in Eq. (H) are qualitatively similar to Eq. ( [4^ 
and we avoid presenting them here. We show the full 
explicit expression for all rotational Love numbers in a 
Mathematica® notebook presented as Supplemental 
Material to this paper. 


EoS Reference 


APR 


MSI 

m 

EPS 

[65] 

SLy4 

m 

BCSK20 


POLYnl P{p) ~ 


TABLE I. List of tabulated EoS used in this work. 


we discuss in Appendix]^ a quadrupolar magnetic tidal 
field would induce a nonvanishing mass moment M 3 at 
linear order in the spin, whereas an octupolar electric 
tidal field would induce a current moment S 2 - 

Furthermore, also the standard moments associated 
with the equatorial symmetry would acquire tidal con¬ 
tributions proportional to the spin. For example, the 
mass quadrupole moment M 2 would acquire 0(x) cor¬ 
rections proportional to the £ = 3 magnetic component 
of the external tidal field (and would also acquire O(x^) 
corrections proportional to the £ = 4 electric component 
of the external tidal field, cf. Paper I). 

An interesting question is whether such deformation 
are still “universal”, in the sense that the relation among 
different multipole moments is still only mildly depen¬ 


dent on the EoS. This problem is discussed in Sec. |V C 
below. 


V. NUMERICAL RESULTS 


C. Deformed three-hair relation for tidally 
distorted NSs 


As anticipated in Paper I, the corrections to the mul¬ 
tipole moments of a tidally distorted spinning NS will 
modify the approximate three-hair relations that exist 
for isolated NSs [59H62]. For such stars these relations 
read [50] 


Mi + i^Si = , (43) 


where a = S'l/M, iq = y^M^TM, Bni is a coefficient 
that depends on the EoS only mildly, and [x] denotes 
the largest integer not exceeding x. Because of the mild 
dependence on the EoS, all moments with £ > 2 can be 
approximately computed from Mq = M, Si = J and M 2 
through relation (43) 


As a result of an external tidal field, the relation (431 is 
deformed in two ways. First of all, we note that such re¬ 
lation is valid for axisymmetric and reflection symmetric 
solutions. An electric (respectively, magnetic) tidal field 
with odd (respectively, even) values of £ would break the 
equatorial symmetry, introducing moments such as 82 ^ 
M 3 , Si, M 3 , etc, which are not included in Eq. (43). As 


In this section, we solve the perturbation equations ex¬ 
plicitly for a self-gravitating perfect-fluid configuration. 
We use some realistic tabulated EoS listed in Table U 
which cover a wide range of NS deformability. We also 
consider a polytropic equation of state, P{p) ~ pi+i/” 
with index n = 1, which will be denoted by “POLYnl” 
in the following. 

Figure shows the mass-radius relation and the mo¬ 
ment of inertia as a function of the NS mass for the un¬ 
perturbed models. As usual, the stellar mass is defined 
by the value of the radial function Ai at the radius, i.e. 
M = M.{R), where the radius R is defined by the rela¬ 
tion P{R) = 0. The moment of inertia reads I = J/0, 
where we recall that is the fluid angular velocity and 
J is the angular momentum. To linear order in the spin, 
the moment of inertia is a 0 {x^) quantity which depends 
only on the compactness and on the EoS. 


A. Tidally-deformed, nonrotating stars 

We now consider the tidal deformations of the back¬ 
ground solution just presented. Let us start by discussing 
nonrotating configurations, which have been studied in 
the past by several authors (e.g. Refs [1 Ellin]) and can 
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FIG. 3. (color online). Left panel: mass-radius relation for the slowly-spinning, perfect-fluid stellar models analyzed in this 
work. We consider various tabulated EoS covering a wide range of NS deformability (cf. Table [I|. Right panel: moment of 
inertia as a function of the stellar mass. 


therefore be used to test our numerical code. 

We solved the perturbation equations for £ = 2, 3,4 
so that we can evaluate the tidal Love numbers in 
Eq. ([^ for £ = 2,3,4. These are shown in Fig. 
where the left and right panels refer to electric and 
magnetic tidal Love numbers, respectively. For com¬ 
pleteness, we show both the tidal Love numbers in 
Eq. (34) (namely A^^) which are associated 

with equatorial-symmetric tidal perturbations and those 
defined in Eq. (E3| (namely, A^\ A^\ A^^) which in¬ 
stead break this symmetry. As previously mentioned, 
these numbers have been normalized by suitable pow¬ 
ers of M to make them dimensionless. Furthermore, in 
Fig. a we multiply the dimensionless Love numbers by 
some suitable power of the compactness C = M/R so 
that, in the Newtonian limit, these quantities tend to 
a constant value. The small-compactness limit can be 
better appreciated for the POLYnl curves, because the 
realistic tabulated EoS are meant to describe only con¬ 
figurations with large compactness. This is the reason 
why the axes of Fig. ^ have been truncated at C ^ 0.05. 
Were these plots extended to C —?► 0, the electric Love 
numbers for POLYnl would approach their Newtonian 
constant value, whereas the magnetic Love numbers all 
go to zero in the Newtonian limit, since axial perturba¬ 
tions are not defined in that case. In the static case, our 
results for the electric and magnetic tidal Love numbers 
are in perfect agreement with those computed in Ref. [6] , 
cf. Appendix [F| for a comment about the computation of 
the magnetic Love numbers in the static case. 


B. Rotational tidal Love numbers 

The rotational tidal Love numbers in Eq. ([^ for var¬ 
ious families of spinning NSs are shown in Fig. which 
is one of the main results of this work. As a reference, 
in Table |TT] we present some rotational tidal Love num¬ 
bers for selected values of the compactness and for two 


c 

APS EoS 

0.1 0.15 0.2 

MSI EoS 

0.1 0.15 0.2 

5X^m^C^ 

5.1508 2.5522 1.2525 

0.0016 0.0024 0.0026 

0.4121 0.1982 0.0930 

0.0046 0.0098 0.0099 

5.9280 3.0727 1.5413 

0.0017 0.0027 0.0031 

0.4789 0.2444 0.1163 

0.0095 0.0160 0.0143 


TABLE II. Selected (normalized) rotational Love numbers 
shown in Fig. for two EoS (APS and MSI) and selected 
values of the compactness. 


representative EoS. Here we only present the rotational 
Love numbers associated with equatorial-symmetric per¬ 
turbations; the case in which the equatorial symmetry is 
broken in discussed in Appendix (see Fig. |^. 

In Fig.[^ we have normalized the rotational Love num¬ 
bers by suitable powers of the compactness C in order for 
the Newtonian limit to be independent of C. The scaling 
has been achieved by investigating the polytropic EoS 
POLYnl as C —>■ 0, since the other tabulated EoS are 
not realistic in the Newtonian limit. 

It is important to note that i5A^^^ ^ 1/C® in the New¬ 
tonian limit, precisely as the static tidal Love number 
A^^ For this reason, as already shown in Fig. ^ the 

ratio JAg does not depend on C for small com¬ 

pactness. 

The left (respectively, right) panels of Fig. (and Fig.[^ 
in Appendix!^ show the “magnetic-led” (respectively, 
“electric-led” Jrotational tidal Love numbers correspond¬ 
ing to a tidal field with harmonic index £ -|-1 (top panels) 
and £ — 1 (bottom panels) that sources a £-pole moment 
with opposite parity, in agreement with the spin-selection 
rules previously discussed. 

Finally, it is interesting to note that not only is 

nonmonotonic, but also it does not have a definite sign as 
a function of the compactness. This peculiar behavior is 
observed also for the other rotational Love numbers as¬ 
sociated with the broken equatorial symmetry discussed 
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M/R 


M/R 


FIG. 4. (color online). Dimensionless electric (left panels) and magnetic (right panels) Love numbers for ^ = 2, 3, 4 (from top 
to bottom) for various nonrotating NS models with different EoS and as a function of the stellar compactness C = M/R. Note 
that, according to our definitions, the electric Love numbers and the magnetic Love number X\,/ are negative. Nonetheless, 
modulo a negative factor, our definitions are fully equivalent to the usual ones (cf. e.g. Ref. [^). Note that A^\ 
associated with equatorial-symmetric tidal perturbations, whereas A^\ break this symmetry. 


in Appendix [cf. Fig. [^ . 


C. Approximate universality 

Some years ago, Yagi and Yunes m HZ] discovered 
that certain suitably normalized combinations of the 
tidal Love numbers of a nonspinning NS and of the multi¬ 
pole moments of a spinning NS are related to each other 
by approximately universal relations which depend only 
mildly on the EoS. In their simplest version, these rela¬ 
tions are known as I-have-Q and connect the moment 
of inertia I, the (rotation-induced) mass quadrupole mo¬ 
ment Q, and the electric quadrupolar tidal Love number 
(all of them suitably normalized) through relations which 
are independent of the NS EoS at the percent level. Sim¬ 
ilar relations also exist among other (electric and mag¬ 
netic) tidal Love numbers of a nonspinning NS jZHj) al¬ 
though the universality deteriorates for increasing values 
of the multipole i and in the magnetic sector relative to 
the electric one with the same i. 

The original I-Love-Q relations were found for slowly, 
rigidly rotating, barotropic, isotropic, unmagnetized and 
isolated stars HSl HZj- Recently, they have been ex¬ 
tended to include rapid rotation (55] [511 [SHI IMj) non- 
barotropic [ZQ|, anisotropic m, and differentially rotat¬ 
ing m fluids, strong magnetic fields m, dynamical con¬ 
figurations [20], exotic compact objects m. and even 
deviations from general relativity [tni HU 1751175] . The 


outcome of these studies is that the approximate univer¬ 
sality is remarkably robust in realistic configurations. 

It is important to stress that the tidal Love numbers 
entering the J-Love-Q relations are those of nonspin¬ 
ning NSs, whereas the moment of inertia and the mass 
quadrupole moment refer to rotating stellar configura¬ 
tions. It is therefore important to understand how spin 
corrections to the tidal Love numbers impact such uni¬ 
versality. The spin corrections to the tidal Love numbers 
of a NS are computed in this work for the first time. 


1. Nonrotating case 

As a by-product of our analysis, we can first study 
the universality of the relations among the (dimension¬ 
less) tidal Love numbers to zeroth order in the spin. The 
results shown in Figs. and confirm the analysis of 
Ref. [55] (but see Appendix|^below for a comment on the 
computation of the magnetic Love numbers performed in 
Ref. [28]) and extend it by computing and A^^ 

In the three panels of Fig. ]^ we show the relations 
among A^^ and A^\ i.e. between magnetic and electric 
tidal Love numbers with the same multipole £. In this 
case the degree of universality is roughly the same for dif¬ 
ferent values of i, typically below the percent level. How¬ 
ever, as noted in Ref. [28j , when written more usefully in 
terms of the lowest multipole moments (for example in 
terms of the dimensionless moment of inertia I = I/M^) 
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M/R 


M/R 


FIG. 5. (color online). Dimensionless rotational electric (left panels) and magnetic (right panels) tidal Love numbers for a 
spinning NS as a function of the compactness M/R and for the various tabulated EoS adopted in this work. In our notation, 
denotes the correction to arising through the coupling to the magnetic octupolar {£ = 3) tidal field to first order 
in the spin. Note that would represent a tidally induced spin shift that can be reabsorbed in the dehnition of x s'-nd 

therefore 5A^^' = 0. 





FIG. 6. (color online). Universal relations for electric and magnetic Love numbers of a nonspinning NS with ^ = 2, 3,4 (cf. also 
Ref. |28|L For each panel, we show the percentual deviation from a EoS-independent fit. The deviations from universality are 
roughly less than 1% between electric and magnetic Love numbers with same multipolar index i. 


the degree of universality deteriorates with increasing £, 
as shown in Fig. for the electric (top panels) and mag¬ 
netic (bottom panels) tidal Love numbers with £ = 2,3, 4. 
For £ = 3 and £ = 4, the universality of the X^ -I and 
relations holds approximately within 5% and 10%, 
respectively. 


2. Universality-breaking for rotational tidal Love numbers 

We are now in a position to investigate the universality 
of the rotational Love numbers computed in this paper. 


In Fig. we show some rotational Love numbers al¬ 
ready presented in Fig.[^but normalized as explained be¬ 
low Eq. ( [M| ) as functions of the dimensionless moment of 
inertia /. The top and bottom panels show the magnetic- 
led and the electric-led rotational Love numbers, respec¬ 
tively. Although the normalized rotational Love numbers 
are less sensitive to the EoS, the deviations from univer¬ 
sality are much larger than in the nonspinning case. For 
example, the rotational Love numbers <5A^ ' (which is 
related to the deformation of the mass quadrupole mo¬ 
ment) shown in Fig. [^displays deviations at the level of 
20% and as large as 50%, to be compared with the de- 




























14 



FIG. 7. (color online). Same as in Fig. [^but for the universal relations between the dimensionless moment of inertia 1 = 1 /M^ 
and the tidal Love number (top panels) or (bottom panels). With respect to Fig. the deviations from universality 
increase with i. 


viations of shown in Fig. which are at the level of 
1 % or better. 

The electric-led case (bottom panels of Fig.l^ is even 
more dramatic, with SX\f displaying deviations larger 
than 200%. Furthermore, as shown in Fig. in this case 
the Love number does not have a definite sign and the 
zero crossing is also EoS dependent. This poses a serious 
limitation to the construction of a fitting formula, since 
at the zero crossing the relative amplitude of SX^^^ for 
one EoS would be vanishingly small relative to the ampli¬ 
tude of the same quantity for a different EoS. This is the 
reason for the “cusps” shown in the right bottom panel 
of Fig. which signal a zero crossing in a logarithmic 
scale. 

For completeness, in Fig. of Appendix [E| we present 
the same analysis but for the rotational Love numbers 
associated with the breaking of the equatorial symmetry. 
The results are qualitatively similar to those presented 
in this section. 

For simplicity, we only present here the fitting for- 
mula for the most relevant quantity, namely (5A^ , the 
fits for other rotational Love numbers are given in the 
Mathematica® notebook in the Supplemental Mate¬ 
rial. The deviations shown in the right top panel of Fig.[^ 


no = -471.133 
ai = -bl358.55 
aa = -1750.22 
as = -bl327.18 
04 = -651.555 
as = -b216.063 


06 = -49.0000 
07 = -b7.50706 
og = -0.74413 
ag = -b0.04314 
oio = -0.00111 


TABLE III. Coefficients of the fit (44 1 . 


are relative to the following fit: 

10 

log <5^ Oi [log /] *, (44) 

i=0 

where the constants Ui are given in Table [lll| As shown in 
Fig. § for the EoS considered in this work (which cover a 
wide range of NS deformability) the deviations from the 
fit ( [44| ) are typically below 30%. 

Although the relation between the rotational Love 
numbers and the moment of inertia depends more 
strongly on the EoS than in the case of the static Love 
numbers, the quantities presented in Fig.j^are first-order 
corrections in the spin, i.e., they are subdominant relative 
to the static case. It is therefore interesting to quantify 
how this variability affects the approximate universality 
of the multipole moments of a NS. To this end, let us 
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FIG. 8. (color online). Spin-induced electric and magnetic Love numbers for spinning NSs as functions of the dimensionless 
moment of inertia 7, to be compared to the static case shown in Fig. Note that the deviations from universality are much 
larger than in the static case. 


consider Eq. (39), which expresses the deformed mass 
quadrupole moment in the presence of an electric 7 = 2 
component and a magnetic 7 = 3 component of the tidal 
field. By defining (A) as the relative spread of a quantity 
A given by different EoS, from Eq. (39) we obtain 




1 % 



where in the last step we used Eqs. ([^ and (1^, the es¬ 
timates (A^) ^ 1%, ^ 40% and JA^^^ 30A^^ 

[cf. Fig. The values for y, ro and C adopted in the 
equation above refer to a typical NS-NS merger. The 
contribution of 5A^ is larger as the binary approaches 


the merger and the approximate universality of the in¬ 
duced mass quadrupole moment deteriorates from 1% to 
roughly 6% due to spin corrections. Such corrections 
scale as x(7?/ro)^/^, so they would be more relevant for 
fastly spinning NSs and negligible at a large orbital dis¬ 
tance. 


VI. CONCLUSIONS AND PROSPECTS 


By extending the formalism developed in Paper I we 
have computed, for the first time to our knowledge, the 
tidal Love numbers of a spinning NS to first order in 
the angular moment. The spin of the object introduces 
couplings between electric and magnetic distortion and 
new classes of rotational tidal Love numbers emerge. 
We have derived the equations governing electric-led and 
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magnetic-led tidal perturbations for generic multipoles 
and have explicitly solved them for various tabulated EoS 
[cf. Table |l] for a tidal field with both an electric and a 
magnetic component with £ = 2,3, 4. 

We found various interesting features; for instance, 
some rotational Love numbers are nonmonotonic func¬ 
tions of the compactness and do not have a definite sign. 
For a binary system close to the merger, various com¬ 
ponents of the tidal field become relevant. In this case 
we find that an octupolar magnetic tidal field can signifi¬ 
cantly modify the mass quadrupole moment of a neutron 
star. Preliminary estimates, assuming a spin parameter 
X ~ 0.05 and an orbital distance tq « 5R, show modi¬ 
fications > 10% relative to the static case. Our results 
suggest that spin-tidal couplings can introduce important 
corrections to the gravitational waveforms of a spinning 
NS binary system. 

Furthermore, some of these new rotational Love num¬ 
bers deviate from a nearly universal relation |26H28) by 
an amount as large as 200% in the electric-led case and 
as large as 50% in the magnetic-led case. By using these 
results, we estimated that - for a NS binary approaching 
the merger - the approximate universality of the induced 
mass quadrupole moment deteriorates from 1% in the 
static case to roughly 6% when y « 0.05 and rg « 5R. 
Clearly, the effect of the spin would be more relevant for 
NSs with X ^ 0.05 and it can only be relevant in the final 
phase of the inspiral. 

The phenomenology of the tidal deformations of a spin¬ 
ning object is very rich and much work remains to be 
done. From the technical side, it is natural to extend our 
study to the nonaxisymmetric case [5], which would in¬ 
troduce a further class of rotational Love numbers emerg¬ 
ing from novel couplings between nonaxisymmetric per¬ 
turbations [cf. diagram in Fig. [^. Likewise, it would 
be interesting (albeit technically challenging) to include 
quadratic corrections in the NS angular momentum, thus 
also extending our analysis in Paper I to the NS case. 

From the phenomenological side, it is important to 
understand spin-tidal effects for the gravitational wave¬ 
forms. This is particularly urgent given that NS binaries 
are the main target of advanced gravitational-wave detec¬ 
tors. Our results suggest that spin-tidal effects might be 
important to improve gravitational-wave templates |79L 
m and to estimate precisely the tidal Love numbers 
through gravitational-wave detections [SHU m Hi]. 

We have shown that the main spin-tidal contribution 
to first order in the spin and in the axisymmetric case 
comes from a deformation of the mass quadrupole mo¬ 
ment of the NS induced by an external octupolar tidal 
field. For a binary system, such an effect is negligible at 
a large orbital distance but becomes more relevant close 
to the merger. Therefore, a solid qualitative outcome of 
our analysis is that spin-tidal effects are more important 
for high-frequency signals, possibly those in which the 
merger phase happens within the detector band. 

A promising approach to incorporate such effects in 
the relativistic two-body dynamics is by extending effec¬ 


tive point-particle techniques (cf. e.g. Refs. to 

include spin-tidal couplings. Work in this direction is 
underway. 

Finally, in Paper I we discussed some subtleties in the 
definition of the tidal Love number for a relativistic com¬ 
pact object (see also Refs. jinilHS]). In brief, in a non¬ 
linear theory the very notion of the multipole moments 
of a single deformed compact object is problematic. We 
believe that a more rigorous analysis is needed to clar¬ 
ify this important issue, even in the static case. This 
would likely require a (at least partial) fifth order post- 
Newtonian expansion of the field equations for a binary 
system, which is currently not available for generic mass 
ratios, in order to determine the quantities which actually 
appear in the post-Newtonian gravitational waveforms. 
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Appendix A: Stationary perturbations of a slowly 
rotating relativistic object 


As a background, we consider the spinning geome¬ 
try © to first order in the rotation rate, with fluid ve¬ 
locity = e“^/^(l, 0, 0, Q) and stress-energy tensor 
given in Eq. ( |Io| ) (see Sec. Ill Aj ). We perform a harmonic 
decomposition of the (stationary) metric perturbations of 
this background in the Regge-Wheeler gauge [86] as 


%.(t, r, £}, g,) = Sg;t^{r, d, ip) + d, ^) (Al) 

with 




/o 0 

* 0 

* * 0 

* * 


0 

0 / 


(A2) 





0 

0 ^ 



0 

0 


^ * 


0 


* 

sin^ J 

(A3) 
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where asterisks represent symmetric components, = 
—1 —2Al/r, = Y^{'d,ip) are the scalar spher¬ 

ical harmonics and we have defined the vector spherical 
harmonics as 


order in the perturbations and in the spin. As discussed 
in the main text, we restrict to a static fluid and therefore 
impose Su^ = 0 in the above decomposition. 


{SlSl,)=(-^smdYU. (A4) 

^ \ sm w ’I 

Here and in the following, a sum over the harmonic in¬ 
dices £ > 1 and m (such that \m\ < £) is implicit. The 
perturbation of the stress-energy tensor reads [40] 

= {P + p){u^6v!^ + + Sp^^'>)Y^u^^u'' 

+6gf^^P + 5P^^'>Y^gi ^'', (A5) 

where 5P^^^ and 5p^^^ are radial functions. The pertur¬ 
bation of the fluid velocity is 


_.(i(pWs^ + i/Wr^^) 

Su = e 2 , 

k{P + p) 

(A6) 

6 u =62 —--- , 

K(P-bp) 

(A7) 

® Atr2(P + p) ’ 

(AS) 

... + 

Atr^ sin^ i?(P -b p) 

(A9) 


Appendix B: Coefficients and sources for the 
inhomogeneous equations governing tidal 
deformations 


The coefficients of the {£ > 2) differential operators 
appearing in the systems and @ read 


(f) _ r {Kr^{P - p) +2)-2M 

Or - 


r(r — 2Af) ’ 

= [r'^cl{r — 2Af)^] ^ [—£{£ + ^)rcl{r — 2M) 
-cl {2Kr^MiUP + bp) + 4M^ 

+Kr"^ (4Kr^P^ — 9P — bp)) 

+Kr^(r — 2M)(P + p)] , 


(Bl) 


_ Kr^(^ + p) 

“ r-2M 


(~<*W _ 

° r(r — 2 Ai 


1 


£{£+!)- — 
r 


Y2Kr‘^{P + p)] , 


(B2) 

(B3) 

(B4) 


where U'^^\ are purely radial functions. The where we recall that Cg = \fdPjWp is the speed of sound 

function has been fixed to ensure = — 1 to first in the fluid. The source terms read 


S' 


* (^) 
-I- 


[l{l -b 1) - 2)^4t{l + 2) + 3c2(r - 2M)2 L 

-cl (2rM (-Kr^p ((b£(£ + 1) - 22) fl + (£(b£ -b 9) - 6)a;i) 


2M)(P + p)(u;i+n)ff^^^ 

-2(£ (£(£ + 1) - 4) - 2) (H-wi) 


+Kr^P (H (3£(£ -b 5) -b 32Kr^p + 14) -b (14 - £(13£ + 25))wi) -b 24K^r^nP^ + SK^r^fip^) 

+r^ (Kr^n (p (b£(£ -b 1) + IGkc^P (kt^P - l) - U) + P (4Kr^P (£(£ -b 3) -b 4Kr^P - 2) + (£ - 3)£ - 6 ) 
-SAcr^p^) + (P (-4k£(£ + 3)r^P -b £{9£ + 13) - 14) -b {£{b£ -b 9) - 6 )p) - 2£{£{£ -b 1) - 2)) 

+2£{£{t -b 1) - 2)H) -b 4M^ {{{£ - 1)£ - 4)(H - uji) - 4Kr^il{P + p))) 


+2r{r - 2M)H^^'’' {{£(£ -b 1) - 2)r(wi - H) + M {{£{3£ -b 5) - 4)(H - uji) + 4Kr^n{P + p)) 

+Kr^P {£{£ + 3){fl — uji) + 4Kr^np) + 4K^r^£lP^))] , (B5) 
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5*^^) = - , ^ - [c^ (2r{r-2M)Hlf'>' (M ((3£(f+ l)-6) - uji) + + p)) 

{£{£ +1) - 2)V4£2 - lc2(r - 2 My L « V ^ ) o \ \\ \ J ) \ 

+r (Acr^P {{^{l - 1) - 2) (11 - wi) + - (£(£ + 1) - 2)(11 - wi) + 4/tV^llP^)) 

( 47^2 ((£2 3 ^ _ 2 ) _ 4Kr2llP - 4«;r2llp) + 2rM {nr^P (- (13£(^ + 1) - 26) wi 

+ (3^2 - 9£ + 2) H + S2Kr‘^np) - nr'^p ((5£(f + 1) - 10) wi + + 9£ - 18) H) 

+2 (^3 + 2^2 _ 3£ _ 2) (H - wi) + 24K2r4llp2 ^2 (-^^2^ ^g^2 

+ (£2 5£ - 2) H - UKr^Qp) + 4K2r^p2 + i) + 2 ) wi + {£{£ - 1) - 4) H + Anr'^np) 

+Kr‘^p {{5i{e + 1) - 10) wi + (5£^ +91- 10 ) n)-2 (f + 2l(i - 1) - 2) (H - wi) + len^r^QP^ - SK^r'^+tp^))) 
+K{£{£+l)-2y{r-2M){P +p){uJi+Vl)H^^A , (B6) 

= . - \(?g (r{r — 2A4)/in^^ (—2A4 (ruji {£(£ + 1) + 2 Kr‘^P + 2Kr^p + 3) 

+ V4FTW+3r2c2(r-2A4)2 L n ^ ^ V iVV A F ) 

+uji (3£{£ + 1) — Anr^P — Anr'^p') — ?>£{£ + 1)11) + r {2K.r‘^P (^£(£ + 1)11 — oji {£{£ + 1) — Anr^p')') + £‘^ruji' 

—2£^n + £ruji + 2£(£ + l)wi — 2£H — 2 F?r'^P'^ (rwi^ — Auji) + 2Kr^puJi) + lOAl^wi') 

+/ig^^ y (^nr'^P (2u)i {£{£ + 1) + l2Kr‘^p) — 2uji' (^r — Anr^p + r) — 5£{£ + 1)H) + Kr^p (2£(£ + l)uji 

—h£(£ + 1)H — Aruji) + 2£{£ + 1)^(0 — wi) — 20K^r®P^u;i — 4K^r^P^ (wi (bnr'^p — A) — Zruji) + 8K^r^p^wi) 

-2rM {-2£^uji + 2£^+l + nr^P (2wi {£{£ + 1) + 22 Kr‘^p + 2) + 2ruji' {Anr^p - 5 - £) - U(£ + l)ll) 

+Kr^p (2 {£{£ + 1) + 2) wi - h£(£ + 1)11 - 12™/) - 21 ^uj^ + 21^VL + £ruji' + Un^r'^P^ (rw/ + 3wi) 

+SK^r^p^uJi — Srwi') — 4A4^ (rwi' (—£ + Snr'^P + Skv^p + lO) + uji {—2£{l + 1) + nr^P + nr^p) + 2£{£ + 1)H) 

+40A4^wi')) — Kr^{P + p)hQ^ (—2rA4 (£{£ + 1)11 — Anr^Pcoi) + (l(l + 1)H + Ak^t'^P^uji) + AM'^uji) , (B7) 

S^- = ^^^2 _ 1 X 2 (r - 2A4)2 {-Kr^P (2wi {£{£ + 1) + Unr'^p) + 2ru}i {£ + Anr^p) - U{£ + l)ll) 

+2£'^{£ + l){n - wi) + Kr^p (-21(1 + l)wi + 51(1 + 1)11 + AruJi) + 20K^r^P^uJi 

PAk^t'^P^ (wi (bur^p — 4) — Srtoi) — S^^r^p^Wi) + 2rA4 (21^a;i — 21^11 + nr^P (2a;i (1(1 + 1) + 22/tr^p + 2) 
+2rwi' (1 + Anr'^p -A)- 51(1 + 1)12) + nr'^p (2 (1(1 + 1) + 2) wi - 51(1 + 1)11 - 12rwi') + 4l2(wi - H) 

—Irwi' + 21wi — 2111 + 12K^r^P^ {ruJi + 3u;i) + — 6rwi') + 4A4^ (ruJi (l + Snr^P + Snr^p + ll) 

+uji (—21(1 + 1) + nr^P + nr^p) + 21(1 + l)ll) — 40A4^wi') — r(r — 2A4)/ig^^ (—2A4 (rwi' (1(1 + 1) 

+2Kr^(P + p) + 3 ) + wi (31(1 + 1) — Anr^P — Anr'^p) — 31(1 + l)ll) 

+r (2Kr^P (l(l + 1)11 — wi (l(l + 1) — Akt'^p)) + £^roji — 21^11 + £ru>i + 21(1 + l)uji — 2111 
—2K^r^P^ {ruJi — Aoji) + 2K.r^ poji) + lOAl^w))) + Kr'^{P + p)^'^ (—2rA4 (1(1 + 1)11 — Anr^Puji) 

+r^ (1(1 + 1)11 + Ak^t^P^uji) + AM^uJi)] . (B8) 

I 

1. Other metric components of the polar sector can be algebraically written in terms 

In their reduced form, the electric-led and magnetic-led 
systems only depend on the variables Hq \ (5/ig^^^^ and 
h\P, 5P[y'^\ respectively. Once such quantities are com¬ 
puted by solving the corresponding system of ODEs, all 
other nonvanishing metric components can be obtained 
through algebraic relations. For completeness we give 
such relations here. 


To zero order in the spin, the nonvanishing components 
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of the function Hq and its derivatives: 




H, 




[{2rM i-L{i - 1) 


(£(e + 1) - 2)r(r - 2M) 

+ 6 Kr^P + 2Kr^p + 4) + {l{t + 1) + 

-2Kr^P - 2Kr‘^p - 2 ) - AM^) 


+2r{r - 2M) (M + nr^P) 

(^) _ Tj(^) 


LJK^) _ ZT 
-^2 ~ -^0 ’ 


dpW = i(p + p)i/W 


(B9) 

(BIO) 

(Bll) 


(i) 

and we stress that i?) = 0 as a results of the field equa¬ 

tions in the stationary case. Likewise, = 0, so that 
the only axial variable in the stationary case to zeroth 
order in the spin is 


To first order in the spin, the nonvanishing pertur¬ 
bation variables which are algebraically related to the 
dynamical ones are 


Sk(^+i) = _ 


nr 


Me 


£{£ + i) [KVWTW+3r^ I 

- 2 VdFTsTTSe'^ diLo ^ ^' 


—-b £{£ + l)uiihQ^ 


2h!i'^' {rui {£{£ -b 1 ) -b Akt'^P) -b 2£{l -b l)a;i - 2l(£ -b l)n) 

2 KrP (2£{£ + 1){£1 - -b V^£‘^ + 8 £ + 

AM^e-‘'uji'h’'^^' , A£{£+l)£le-‘'iP + p)k 


kVA£^ -b 8^ -b 3 L 


-b 


^0 




KV4FTWT3r2 V4£2 -b -b 3 

[—2A4^ (—4a;i (£(£ + 1) — nr'^P — nr'^p) + roji {—£{£ + 1) + Skv^P — 2) 


kV 4£2 -b8£-b3r3(r - 2M) 

+A£{£ + 1)12) - (wi {£‘^ {£^ -b 4£ -b 3) -b ^K^r^P^p -b Sk^t^P^) - {£^ -b 4£ -b 3) 

-bwi' (k£{£ + S)r^P -b 2£r — Ak^t^P^)) + rM (roJi' {2k -b 3^ -b 4) r^P — {£— £>)£ — Sk^t'^P^) 
-b2wi {£ {£^ + A£{£ -b 1) - 2) - ^K^r'^Pp - %K^r^P^) - 2£ (£^ + A£{£ -b 1) - 2) 12) - 

(2rM (-£^ - 3£ + Akt^P + 2)+r^ {£{£ -b 3) -b Ak^v^P^) - 4A42) ^ ^ ^ ^ ^ ^ 

Kr^{r-2M) ^ ^ 


(B12) 


= _ 


{£{£ - 1) - 2) V4F^r(r - 2M) L 


-r(r - 2M) (2M (4^^' (rw/ {£{£ -b 1) -b Akt^P) + 2 £{£ + l)(wi - 12)) 


-b\/4£2 - (^-2KrP (2£{£ + 1)(12 - wi)4^^' - -b 4 kV^P2^i'4^^' 

-£{£ + l)wi'/i[,^^') -b AM^uji%^o^'^ + 2hi^'> (2M^ (Auji (£(£ + 1) — Kr'^P — KT^ p) + ruji {£^ P3£ — Skt^P -b 4) 
-A£{£ + 1)12) - rM {-2Kr^P {{£(£ - 1) -b 2) rw/ -b 2£{£ + 1)12 - S^r^pwi) 

+{£ + 1) (2f3l2 -b 2 (4 — £(£ -b 1)) £oJi — 2£‘^Q — AK£r^£lp + £ru)i — 8^12 -b 6rwi') -b SK^r^P"^ {rwi -b 2wi)) 

-r^ {-{£ -b 1) (£312 -b {-£{£ -b 1) -b 2) £wi - £'^n - 2K£r^l2p - 2£12 -b 2ra;i') 

-bK(£ -b l)r^P ((£ — 2)ruji -b 2£12) -b Sn^r^P^coi + An^r^P'^ {2Krpuji — wi')) — 8Al^wi') 

-e^V4£2 - (2rM (-£(£ -b 1) -b 6Kr^P + 2Kr^p -b 4) -b (£(£ - 1) 

PAk^t^P^ - 2Kr^(P + p)-2) - AM"^)] , (B13) 
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5H. 


(^+1) 


(-W 1 (£(€ + 1) + 4.Kr^{P + p)) + £{£ + l)f2 + 2(r - 2M)uJi') - 2r(r - 2 M)uJi'h^o^ 

(B14) 


^M{£ + 2) + 3 L' 

+ ^4i{i + 2) + 3e''SHo^‘^^^^ 

= - „ [2^4^^ (-wi (€(£ + 1) + 4Kr^{P + p)) + i{£ + 1)0 + 2 {r - 2M)wi') - 2 £r{r - 2M)a;i'4^^ 


£V4£2 - 1 L' 




(f-i) 




e-^(P + p) 
2y/4£{£ + 2) + 3 L' 


2 £(^+ 1 ) 0 ^^^ + ^/4£{£ + 2 ) + 3e’'SHo’'^^^'> 


= -e-‘'{P + p) i - 


2 £{£+l)nh''Q^' 
\/4£2 - 1 “ 


(B15) 

(B16) 

(B17) 


whereas = 0 = so that the perturbed the main text, such perturbations with £ = 1 generically 

metric is symmetric under t ^ —t and ip ^ —ip. satisfy a different set of equations. 

In the electric-led sector, it turns out that the equa¬ 
tions for £ = 1 can simply be obtained from those pre¬ 
sented above after setting £ = 1. However, the £ = 1 
magnetic-led sector is different. In this case the equa- 
To first order in the spin, a quadrupolar tidal field tions for Pl^'^ have the same schematic form as (191, but 
will source dipolar metric perturbations. As discussed in with 


2. Perturbation equations with £ = 1 


cW = o, 

^(1) _ _ ^ _ 

° c2(r — 2AI)2 (AI + Acr^P) 


[cl (M (kv^ (lOp - Kr‘^{P - p)(19P -H 3p)) - 6) - nrM^iSbP -1- lip) 


+Kr^ {nr'^ (Pp (S/ir^P + 2) - Sp^ (kt'^P + l) + P^ (3 - 2Kr'^P)) + 2p)) 
-|-Kr(P -I- p) (-|-Kr^P + r — M) [M + Kr^P)] , 


(B18) 


(B19) 




5^ 


r^clir - 2M)2 (M ~ 

Pruji — \2{ijJi — O)) -I- [nr^P (rwi' -|- 28wi — 280) -|- p (—Srwi' — dwi -|- 40) — Iruji — 8wi -I- 80)^ 

-|-r^AI {2K^r^P^ ( 2 /ig^^ {rui + bwi -I- O) -|- r/ig^^ {tloi — 12(a;i — O))^ — nrP (2Kr^p {3ru!i + 2uJi — 120) 
+7ruii' - 28wi -h 240) -h {2Kr‘^p (Srw/ - 4uji + 40) - 9rwi' - 28a;i -f 280)) -t 12K^r^Vlp'^h^Q^ 

—nrp (18ra;i' -|- 32a;i — 400) -|- (9rwi' -I- 4a;i — 40)) -|- uji' + lOhg^^)) 

{ 2 h^^\n - uji) - {ruJi' + 2 uji - 20)) -t (leK^r’^OP^hl,^) 

+KrP (nrp (3ra;i' — 4wi -|- 40) — Auji) — 2h[^^ [Gn^r'^Hp^ + ur'^p (—3ra;i' -I- 4wi — 40) -|- 8a;i — 120)) 

+K^r^P^ (2Kr^0p — ruJi + 6O) — rh!'^'’ {ruJi — 12{uJi — O))) 

— (3«:r^p—2) ^l/trOphg^^ — Wi'^rhg^^ -|-2/ig^^)))) -|-4H;r^0hQ^^(P-|-p) (rAI — AI^ + Kr^P («:r^P-|-l)) .(B20) 


When £ = 1, we can use a residual gauge freedom to set 
if(i) = 0 in the ansatz for the metric [mill], so that the 


r 


only further nonvanishing polar components read 


.rrW ^ e 

^ 5 {M + nr^P) 


A4 (2-\/l5r^wi'hg 


{ruji' -£ 3a;i - 30) -£ 15e^(5P^^^) 
-l-r^ (^nrp ^4-\/l50/ig^^ — be'^5 h'^^ 


-l-VTs ( 4KrPuJih[?^ — oji (r/ig^^ -1- 2 h, 


,( 2 )' 


,(2) 
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Appendix C: Exterior solution for ^ = 3 and £ = 4 

The explicit solution describing the exterior of a slowly 
spinning object immersed in an axisymmetric £ = 3 elec¬ 
tric tidal field reads 


^( 3 ) _ 2(2/ - 2)2(y - l)y^a 3 + (2 (-2 - lOy + 65y^ - GO//^ -f lby‘^) -f lh{y - 2f{y - l)y‘^ log [1 - 2/y\) 73 

° “ 2(2/-2)2/ 


(Cl) 


6 h 


(4) 

0 


Mx 

33602/2 


- 64 ^ 72 / {2-9y + 5y^) as + lOV? (2 (I 6 -f 76y - 762y^ 


2610y3 -f 198 652/"^ - 213152/® -f 60 902/®) 


+3y (-32 -f 1442/ - 80y^ - 5800//® -f 130502/“^ - 9135?/® -f 2030//®) log [1 - 2/y]) 73 

-£ 492 / (2 ( 2 IO 2 /® - 4 - 182 / - 90y^ + 685//® - 7352 /^) -£ 152 /® (SO// - 63//^ -£ 14?/® - 40) log [1 - 2/y]) 743] , (C2) 


Sh 


( 2 ) _ 
0 ~ 


Mx 

" 42002/2 

-£4802/^ 


-8^352/ (16 + 402/ - 240//^ -£ 105?/®) 03 + 5 {Ay + 30Ay‘^ - 774//® -£ 354//® - 8 -£ 2 / (80// 

- 5642 /® -£ 177y'^ - 32) log [1 - 2/y\) 73 -£ 35// (3(y - 2)//® log [1 - 2/y\ - 4 - 4// - 6 //^ -£ 6 //®) 723 )] 

(C3) 


where, as in Sec. lllE y = r/M. Likewise, for an axisymmetric £ = 3 magnetic tidal field, we obtain 


hi®) = - 


= 


M (- 642 /® (8 - IO 2 / -£ 32 / 2 ) -£ 3 (8 -£ 20// -£ GO//^ - 210?/® -£ 90//^ -£ 15//® (8 - 10// -£ 3y^) log [1 - 2/y]) j* 

-512^72/^ (18 - 25y + 8 y^) ag -£ V®7 (-192 - 288?/ -£ 480y^ + 4520//® -£ 46 802/'‘ - 219840//® 


04) 


2688(2/ - 2 ) 2 /® L 

-£404200?/®-24 6 7 502/^-£49 3 502/®-£ 15?/^ (432 -£ 5040y - 186 082/2 -£ 21150//® - 9870//^ -£ 1645?/®) log [1 - 2/y]) 7 ^ 
+672//® (2 (4 + 362 / - 4802/2 + 860//® - 525//^ + 105?/®) + 15(?/ - 2 ) 2?/2 (6 - lAy + 7y^) log [1 - 2/y]) 74 * 3 ] , (C5) 


5ij(2) =-_y-^ 

° 2240(2/ - 2 ) 2 /® L 


-2560352/"^ (38-312/-£62/2) a3-3035 (2 (-16-24?/-722/2-902/®-£31202/^-4005?/®-£1215//® 


+ 152 /® ( 4482 /- 3482 / 2 + 8 l 2 /®- 152 ) log [1-2/y]) jl + ll20y^ {4+8y-18y^+6y^+3{y-2fy^ log [l-2/y]) 72 * 3 ] ■ 


(C 6 ) 


Finally, the solution describing the exterior of a slowly spinning object immersed in an axisymmetric £ = 4 electric 
tidal field reads 


= 28{y-2)y ~ ^4 + 7 (2 (4 + 36y - A80y^ 


1 


860//® - 525//® + 105?/®) 


+ 15 ( 2 / - 2)^2/" (6 - 142 / + V) log [1 - 2 /y]) 74] 


(C7) 


<5hi®) = - 


Mx 


-2560112/ (-16+1682/-2402/2+852/®) ^4+35 flGOlI (-48 - 3442/+82 562/2 + 3 19802/®-45 9 2 702/ 


4435202/2 
+884520//® - 586845//' 

+42525//^) tanh- 


1 


127575//^ + 3// (128 - 1344// + 19202/2 + 112720//® - 396900//® H 
74 + 99y (2 (8 + 56// + 4202/2 - 56 70//® + 10920//® - 7245//® 


476280?/® - 238140//® 


1-2/. 

+ 105?/® (40 - 140// + 1682/2 - 84//® + 15?/®) log [1 - 2/y]) 754 )] , 


1575?/® 


(C8) 


__ Mx _ 

0 28224y2 


-12807// (-4 - 30// + 210?/2 - 224//® + 63y®) ^4 + 7 ^407 (96 + 40// - 123002/2 + 52140//® 


-62370//® + 20250//® + 15?/ (32 + 240// - 16802/2 + 35 92?/® - 2754//® + 675//®) log [1 - 2/y]) 74 + 63?/ (8 + 20// 
+60?/2 - 210//® + 90//® + 15//® (8 - 10// + 32 / 2 ) log [1 - 2/y]) 734)] , (C9) 
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whereas the magnetic counterpart reads 




= - 


M 


[-240y3 (-40 + my - + 49 (2 (-4 - 18y - 90y^ + 685y^ - 735y^ + 210y^) 


= 


3360y 

+15y3 (_40 + goy _ 63j/ 2 + I4y3) log [1 _ 2/2/]) 7 *] , (CIO) 

55440(y - 2 ) 2/4 [-480yTIy4 (-350 + 939y - 732y^ + 175y^) + 7 (^7VTi (2 (24 + 68y - 60y^ - 15622/^ + 6882y^ 

+ 128 8 342/® - 42 1 8902/® + 462735//^ - 213570?/® + 35 5 952/®) + 15/ (-700 - 4450y + 30 1 762/® - 55020//® + 44 2 962/'‘ 
-16611?/® + 2373//®) log [1 - 2 /y]) 74 + 1980//® (-8 - 112?/ + 2576//® - 7560//® + 8190?/^ - 3780//® + 6302/® 

+ 105(2/ - 2)®?/® (-2 + 82 / - 9?/® + 3?/®) log [1 - 2/y]) j*^)] , (Cll) 


X 


480v^ 2/'‘ (-730 + 1401y - 798//® + 140//®) a)] + 49 (/7 (2 (24 + 68 ?/ + 300?/® + 856//® - 36360//^ 

(C12) 


35280(2/ - 2)2/4 L 

+ 102095//® - 80220//® + 19005//^) + 17>y^ (1460 - 7870y + 11732?/® - 6615//® + 12672/4) log [1 - 2/y]) 7 ^ 
-360//® (2 (-2 - lOy + 65?/® - 60//® + 152/4) + 15(2/ - 2)®(2/ - 1)?/® log [1 - 2/y]) 7 * 4 )] . 


The solutions above depend on various integration con¬ 
stants. As explained in the main text, the constants ai 
(respectively, a *) are proportional to the electric (respec¬ 
tively, magnetic) components of the external tidal field, 
whereas the constants 7 ^ and jij (respectively, 7 * and 
7 j* ) are associated with the deformation of the mass (re¬ 
spectively, current) i = i multipole moments. Terms 
proportional to 7 i and 7 * are of zeroth order in the spin, 
whereas terms proportional to 7 ^- and 7 )) are of linear 
order in the spin. 


Appendix D: Tidally-deformed, spinning black-hole 
solution 

Although our interest in this work is devoted to NSs, 
for completeness we present here the explicit solution for 
a tidally-deformed spinning black hole to first order in the 


spin, in the presence of both an electric and a magnetic 
tidal field with £ = 2,3,4. This solution complements 
those presented in Ref. [3] and in Paper I by presenting 
the axisymmetric £ = 3,4 (both electric and magnetic) 
terms to linear order in the spin for the first time. 

As explained in Paper I, when the central object is a 
black hole the requirement of regularity at the horizon 
simplifies the metric perturbations considerably. Using 
the same procedure detailed in Paper I, it is straightfor¬ 
ward to show that regularity at the horizon implies that 
all 7 ’s be exactly zero. Therefore, in agreement with the 
findings of Ref. [2] and of Paper I, all tidal Love numbers 
of a spinning black hole are vanishing. 

The nonvanishing components of the line element 
in the case of a spinning black hole deformed by a station¬ 
ary, axisymmetric tidal field with electric and magnetic 
components with ^ = 2, 3,4 to linear order in the spin 
read 


gtt = ——^ —896v^+ 112 -\/ 5 ( 2 / — 2)y{l + 3 cos 2 ??)a 2 + 56V7y (2 — 3y + 2 /®) (3cos?? + 5 cos 3 ? 9 )q !3 — 324ya4 
896v7r?/ L 

+918//® a 4 — 756//® 04 + 189?/4a4 — 7202/cos 2 ??q ;4 + 20 402/®cos 2 ??a 4 — 16802/®cos 2 ??a 4 + 420?/4cos 2 ??a 4 
— 12602/cos4??a4 + 35 702/®cos4??a4 — 2940//® cos 4 ??a 4 + 73hy‘^cos4'daA 


(y - 2)x 


672/5COS 'd{3-4y + 3cos 2?9)a2 + 28\/7(-55 + 24y + 4(-39 + 16?/)cos 2?? + 5(8// - 9)cos 4'd)a; 


+6cos ^ (5 (-17 - 62?/ + 63?/®) + 4 (675 - 854// + 175//®) cos 2d+ 7 (l75 - 382y + 175?/®) cos 4??) a^] , (Dl) 


9t(p — 


M{y — 2 ) 2 /®sin'd® (336-\/5cos??a 2 + 28-\/7(—4 + 3y){3 + 5 cos 2 ??)a 3 + 15 (20 — 35// + 14//®) (Ocos-d + 7cos3??)a|) 


2 ^ 


224v^ 

112v^ (-12 + 15?/ + 20?/® + 5(5// - 4)cos 2'd) + 56v^ ((22 - 35// - 105//® + 60//®) cos?? 


Mysin d 
2240i/7r2/ L 

+5 (2 - 9?/ + 5//®) cos 3??) a 3 + 5 (896^^ + (9 (-16 + 40// + 240//® - 351//® + 1122 / 4 ) + 4 (-64 + 192?/ + 840?/® 
-1295?/® + 4202/4) cos2?? + 7 (-16 + 168// - 240//® + 85//®) cos4??) a 4 )] , (D2) 
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Qrr — 


896y7r(?/ - 2) L 


896v^ + 112-\/5(y — 2)y(l + 3cos 2 -d)a2 + 56V7y (2 — 3y + y^) (3cosi? + 5 cos 3 '(?)a 3 — 324yQ;4 


+918j/^a4 — 7562 /^a 4 + 189y^a4 — 720?/cos 2 i 5 a 4 + 2040y^cos 2 '(?q ;4 — 1680?/^cos 2 'da4 + 420i/^cos 2 'da4 
— 1260ycos4i?Q;4 + 3570y^cos 4 i 5 a 4 — 2940y^cos 4 i?q ;4 + 735y"^cos 4 z?a 4 ] 


yx 


134475008 i?^a 2 + 2877(-61 + 30y + 4(-45 + 22j/)cos 2^* + 5(-3 + 2?/)cos 4^l)a5 


8967^(2/ - 2) L 

+ 6 cos i? (-5 (147 - 260?/ + 9l7) + 4 (1065 - 1820?/ + 6377) cos 2?? + 21 (15 - 20y + 7?/^) cos 4^ a^] , (D3) 


gm = 


M'^y'^X 
134477? L 


1477(354 - 261y + (1096 - 804?/)cos2?? + 5(-34 + 21 ?/)cos 4 ??)a 3 + 18cos?9 (1369 - 2295y + 8687 


-4 (1175 - 2079?/ + 8127) cos 2?? + 21 (55 - 81y + 28?/^) cos Ad) 04 ] 




56075 (?/2 _ 2 ) (1 + 3cos 2 ??)a 2 + 5677 (4 - lOy^ + 57) (3cos d + 5cos 3 ? 9 )a 3 


44807?? L 

+5 (8967?? + (-8 + 60?/2 - 707 + 2l7) (9 + 20cos 2d + 35cos4??)a4)] , 
= gM sin^ d , 


(D4) 

(D5) 


where we recall that ai and are proportional to the 
amplitude of the electric and magnetic tidal field, respec¬ 
tively [cf. Eq. (30)]. 


Appendix E: Tidal Love numbers with broken 
equatorial symmetry 

In the main text we have focused on tidal perturba¬ 
tions that preserve the reflection symmetry of the back¬ 
ground relative to the equatorial plane at ?? = 7 r/ 2 . A 
binary system of two nonspinning compact objects or of 
two spinning compact objects with (anti)aligned spins 
enjoys such equatorial symmetry relative to the orbital 
plane. However, in other configurations - e.g. if the an¬ 
gular momenta of the two objects are not (anti)aligned - 
the equatorial symmetry can be broken. Supernova kicks 
after the formation of a NS might misalign the spin, 
although subsequent tidal dissipation contributes to re¬ 
align them. The details and the efficiency of these pro¬ 
cesses are largely unknown (cf. Ref. [SS] and references 
therein). For completeness, in this appendix we briefly 
extended the results discussed in the main text also to 
nonequatorial-symmetric tidal perturbations. 

A tidal perturbation is symmetric under reflection 
around the equatorial plane (equatorial-symmetric for 
short) if it is invariant under the transformation d —)■ 
TT — d. Given the symmetry properties of the spherical 
harmonics, electric (polar) perturbations are equatorial- 
symmetric when m is even, whereas magnetic (axial) 
tidal perturbations are equatorial-symmetric when i + m 
is odd. In the axisymmetric case, electric (respectively, 
magnetic) perturbations are equatorial-symmetric when 
i is even (respectively, odd). 

The technique adopted in the main text to extract 
the deformed multipole moments is based on Ryan’s ap¬ 
proach |53H55j and assumes equatorial symmetry. How¬ 


ever, the Geroch-Hansen moments are also defined for 
geometries that break the equatorial symmetry. For such 
geometries the mass multipole moments Mg with odd £, 
and the current multipole moments and Si with even £, 
are different from zero [521 (SO]- Although it should not be 
difficult to extend Ryan’s approach to the nonequatorial- 
symmetric case, this would bring us too far. Therefore, 
here we opt for a more practical definition; a detailed 
analysis will appear elsewhere. 

For our purposes, it is sufficient to notice that the mul¬ 
tipole moments can be extracted (modulo a normaliza¬ 
tion factor) from the large-distance behavior of the metric 
in our coordinates, namely, 

e =2 ^ 

(E2) 

e=2 


In the following we will use the above large-distance be¬ 
havior of gtt and g^p to define the multipole moments 
associated with nonequatorial-symmetric tidal perturba¬ 
tions. Glearly, such a procedure does not allow us to fix 
the overall factors of Mi and Si and it is not manifestly 
gauge invariant. Our procedure here should therefore be 
extended through a more rigorous analysis, but it is nev¬ 
ertheless sufficient to capture the correct scaling of the 
multipole moments with the compactness. 

In the static case, in addition to the Love numbers 


defined in Eq. (|34|, we can also define 

7 ^ A/f — . 

as 


\( 3 ) _ 73 

'^E — „ > 


72 ) = 

^ 480aS’ 


r( 4 ) ^ 74 

^ I1520a7 


(E3) 


where the prefactors have been chosen in order to agree 
with previous definitions of the static electric and mag¬ 
netic Love numbers |1|7|. The Love numbers in Eq. ( |E3| ) 
are shown in Fig. [^in the main text. 
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To linear order in the spin, similar to the derivation of 
Eqs. (35|-(38), we can define the (dimensionless) quan¬ 
tities 




M ^X- 


5\ 


5\ 


74 v ^73 + 35723 


16800a3 


( 32 ) 


oc -X- 


.'M 

a% 


( 34 ) 637^771 - 36073*4 


(X X 


360aJ 


j, (43) 374\/7 73 -I- 49743 

^ 56448003 


(E4) 

(E5) 

(E6) 

(E7) 


which are proportional to the rotational Love numbers 
in Eq. (34) for the nonequatorial-symmetric case. Note 


that a quadrupolar magnetic tidal held (which breaks 
the equatorial symmetry) would source a mass moment 
Ml. However, such dipolar contribution can be elimi¬ 
nated through a shift of the central object along the axis 
of symmetry and is therefore irrelevant. 


10 which can be 


Figs. 1^ and ^^ 

compared to their counterparts (Figs. and m 
the equatorial-symmetric case. We observe that 
the nonequatorial-symmetric rotational Love numbers 
present the same qualitative features discussed for the 
equatorial-symmetric case in the main text. 


The rotational Love numbers dehned in Eqs. (|E4 1- 
(E7) are shown in 
their 


Appendix F: Note on the magnetic Love numbers 
for static fluids 


In this appendix we point out an inconsistency in the 
computation of the magnetic Love numbers for static flu¬ 
ids that was not reported in previous work. Although this 
issue is not directly related to our main goal, here we take 
the opportunity to correct some results that appeared in 
the past. 


In the static case, magnetic tidal perturbations of a 
static fluid are governed by the first equation in (19). 
Using the differential operator (18) whose coefficients are 
given in Appendix]^ this equation explicitly reads 


El 


d'^hQ'^ Kr'^{P + p) 
dj .2 j. _ 2 j\4 dr 


l(i + 1) — 4AI/r -I- 2Kr‘^{P + p) ,{e) 
r(r-2M) 



Such an equation is equivalent to that obtained by Bin- 
nington and Poisson (cf. Eq . (4.29) in Ref. [6]). Although 
Ref. [5] reports that Eq. ( |F1[ ) is also equivalent to Eq. (31) 
in Ref. [7], this is actually not the case. To show this, let 


us first write Eq. (31) in Ref. [7] in our notation: 

2M{r) + Kr^[P(r) — p{r)] 


E. = 


dr'^ r (r — 2Ai) dr 

.2/ 


£(£+l)_6M/r-^r^(P-p) _ 

r{r-2M) ^ U, ( j 

where is a gauge-invariant master function |87) . As 
discussed in Ref. [7], for static perturbations = 

rh^ — 2hg . Clearly, by inserting this definition into 
Eq. (F2) one would obtain a third-order ODE for hg . 
By using the background field equations and Eq. (FI), it 
is straightforward to show that 


El — E 2 = 


Kr{P + p) 
c2(r-2M)2 


[ 2 AK 


(f) 


rBh, 


O' 


(F3) 


where 


A = cl{r {£{£ + 1) -b Kr^{P + 2p) - l) - 3M) 

—M — Kr^P , 

B = cl{r {nr^iP + 2p) + 3) - 7M) - M - nr^P , 

and we recall that Cg is the speed of sound in the fluid. 
The difference Ei — E 2 cannot be further simplified. 
While the two equations are equivalent in vacuum where 
P = 0 = p, they are generically inequivalent inside mat¬ 
ter. In principle, this is not enough to show that the 
magnetic Love numbers defined through these two equa¬ 
tions are also different, since the matching procedure has 
to be performed on different exterior solutions. However, 
as a further proof that Eqs. (El) and (F2| are inequiva¬ 


lent, we have computed the magnetic Love numbers for 
£ = 2 in both cases, by matching the interior solutions 
and with their corresponding exterior solution 
at r = R. In the latter case, we define the magnetic 
Love numbers as in Ref. [7] (see also Ref. [28]) modulo 
an irrelevant prefactor. 

A comparison between the two cases is shown in 
Fig. 11 where the (dimensionless) magnetic Love num¬ 
bers computed from Eq. (F2) have been arbitrarily nor¬ 


malized. As this hgure clearly shows, the dimensionless 
magnetic Love numbers are inequivalent, although nearly 
universal in both cases. We verified that the same result 
holds also for £ = 3 and £ = 4. 

We believe that Eq. (31) in Ref. [7| is incorrect, proba¬ 
bly due to a flawed limiting procedure to the static case. 
As a consequence, recent results which have used that 
equation to compute the magnetic Love numbers (e.g. 
Ref. |2H] , whose qualitative results have anyway been con¬ 
firmed in this work) should be revisited. 

Finally, we remark that recent work by Landry and 
Poisson [49] has shown that the axial perturbation equa¬ 
tions are much more naturally solved in the case of an 
irrotational (rather than a static) fluid, as can be seen by 
taking the static limit of the full perturbation equations. 
The issue just discussed is not related to the corrections 
presented in Ref. |49j . which we hope to include in our 
analysis in the future. 
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FIG. 9. (color online). Same as Fig. but for the case of tidal perturbations that break the equatorial symmetry. Note that 
in this case the rotational Love numbers are defined modulo a prefactor. 
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FIG. 10. (color online). Same as Fig. but for the case of tidal perturbations that break the equatorial symmetry. Note that 
in this case the rotational Love numbers are defined modulo a prefactor. 
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